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Abstract 

The process of electron-positron pair production by linearly polar- 
ized photons is used as a polarimeter to perform mobile measurement 
of linear photon polarization. In the limit of high photon energies, 
w, the distributions of the recoil-electron momentum and azimuthal 
angle do not depend on the photon energy in the laboratory frame. 
We calculate the power corrections of order m/tu to the above dis- 
tributions and estimate the deviation from the asymptotic result for 
various values of uj. 



1 Introduction 

The differential cross-section for electron-positron pair production by linearly 
polarized photons was derived in a series of papers during the period 1970- 
1972 ^ (see also p|, Q and references therein). Expressed as a function of 
s = 2muj, where m is the electron mass (which we shall set equal to unity) 
and 00 is the photon energy in the laboratory reference frame, the differential 
cross-section with respect to the azimuthal angle, 0, between the photon 
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polarization vector e and the plane containing the initial-photon and recoil- 
electron momenta, is given by 



dcf) rn? 

with 



28^ 218 20 
Y ~ ^ ~ U ~ 27 



P = eisin(20)+e3cos(20), i: = ln^. 



Here and ^3 are the Stokes parameters describing the photon polarization, 
introduced through its spin-density matrix: 

Pij = e-e-= i(l + (Ti)ij. 

In the derivation of (|I]) terms of order iv? / s were systematically neglected. 
The main contribution, ~ 0{L), arises from configurations with small re- 
coil momentum, q = fn^^j ^ fn, where 9 is the polar angle of the recoil 
electron (i.e., the angle between the initial photon and recoil-electron direc- 
tions). However, the corresponding events presumably cannot be measured 
experimentally. For the region g ~ m (^^ ~ 50°), the doubly differential 
cross-section was obtained in 0: 

-,9 asym „ o 

with 

. 2£: - 3 , , , , , 1 , / X 

ro = — , ao = IH ln(g + £:), Oo = 1 ln(g + £:), 

m q q 

where e = \/q^~l^. The comparatively large magnitude of the azimuthal 
asymmetry 

A=^ = -- — g2 + _^qi + 0(g6) ^ 14% (3) 

ao 7 245^ 34300^ ^ ; v ; 

is, in fact, the reason this process is used for the polarimetry of linearly 
polarized photons P, 

The aim of the present paper is to calculate the power corrections of 
order 1/s to the asymptotic expression for the asymmetry. The calculation 
of radiative corrections to the asymptotic expression for the asymmetry is a 
rather difficult problem, which we shall not touch here. A rough estimate 
gives A^'"^'^ ^ ^Lr~^2- 3%. 
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The differential cross-section of electron-positron pair photoproduction 
off a free electron in the Born approximation is described by eight Feynman 
diagrams. It was calculated numerically in particular by K. Mork The 
closed expression for the unpolarized case is very cumbersome and was first 
obtained in a complete form by E. Haug during the period 1975-1985 [Q. To 
the best of our knowledge, the exact analytical expression for the differential 
cross-section in the case of a polarized photon has not yet been pubhshed. 
Special attention has been paid to the so-called Bethe-Heitler (BH) subset of 
Feynman diagrams, whose contribution does not vanish in the high-energy 
limit, s ^ oo 1]^. The power corrections to this contribution to the total 
cross-section, behaving as L^/ s [§], indicate the need for the exact expression. 

A detailed analysis of the expressions of Haug's work reveals that the 
interference terms of the BH matrix elements with the other three gauge- 
invariant subsets (which take into account the bremsstrahlung mechanism of 
pair creation and Fermi statistics for fermions) turn out to be of the order of 
some percent for s > 50 — 60m^. On the other hand, the difference between 
the asymptotic and the exact expression is still found to be of the order of 
several percent for s > 3000m^, i.e., very far above threshold, rendering the 
asymptotic expression useless in the energy range of interest. Arguments 
of positivity of the cross-section provide the relevant upper bound for the 
polarized part of the differential cross-section. 

In Ref. a Monte Carlo simulation of the process under consideration 
was performed using the HELAS code, in which all eight lowest order di- 
agrams can be numerically treated without approximation. There it was 
shown that one might consider only the two leading graphs in a wide range 
of photon energies from 50 to 550 MeV. Note that this observation was made 
earlier for the unpolarized case by Haug (who presented his results in 
explicit analytical form). 

Our paper is organized as follows. After introducing the problem, in sec- 
tion 1^ we analyze the kinematics of the process and give a general expression 
for the differential cross-section taking into account only leading and non- 
leading (~ 1/s) contributions. Section ^ is devoted to the derivation of the 
differential cross-section with respect to the azimuthal angle and the recoil 
momentum of the electron. In the concluding section we present the cor- 
rection to the cross-section and asymmetry, together with some numerical 
estimates. Some details of the calculation may be found in the Appendix. 
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Figure 1: The Feynman diagrams contributing to triplet production 

2 Kinematics and differential cross-section 

In the Born approximation, the cross-section for the process of pair produc- 
tion off an electron, 

j{K, e) + e{p) ^ e{q) + e(p_) + e(p+), (4) 

with 

q = K — ki, p^ = p — k, p^ = ki + k, 

is described by eight Feynman diagrams (Fig. |T]), which can be combined into 
four gauge-invariant subsets. Bearing in mind the desired application to the 
case of high photon energies, u; ^ m, we shall present the total differential 
cross-section with leading terms (non-vanishing in the limit s = 2muj 00) 
and terms of order 1/ s (non-leading contributions). The first arise from the 
BH subset, denoted by the indices (12), whereas the non-leading terms come 
from interference of the BH amplitude with the sets denoted (34), (56) and 
(78), as well as from the BH amplitude itself. 

We use the following Sudakov decomposition of the momenta in our prob- 
lem: 

k = akP + PkK + k±, ki = avp + xK k\^, 
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2 

p' =p- K — , s = 2pK, (5) 
s 

with the properties that 

p'2 = = 0, 2p'p = 1, k±p = A;^^ - 0. 
The energy-momentum fractions of the created pair are x and 1 — x while 

— * — * — * 

—ki and k — ki are their momentum components transverse to the photon 
beam axis. The three mass-shell conditions, 

{K-kif-l = -sai{l-x)-A = 0, A = k^ + I, 
{ki + kf - 1 = s{ak + ai){l3k + x) - B = 0, B = {k + k^f + 1, (6) 
{p-kf-l = -spkil - cck) -P-ak = 0, 

permit elimination of the following three Sudakov parameters: 

Pi- A 

Pk = , ak = -{si + k'^), ai = — -. (7) 

s s s{l — x) 

Henceforth, we neglect terms contributing to the cross-section at order ~ 
Keeping in mind the application to experiment, we shall assume 

P ~ 1. (8) 

In this context, it should be noted that the above expressions for ak and f3k 
are also valid only to order ^(s^^). Here Si = {k + KY denotes the invariant 
mass squared of the created pair and /c^, the recoil-electron four- momentum 
transfer squared. 
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{x + Pk){l-x) 



(fcl + (l-x)^)' + l 



k'^ 



The recoil-electron 3-momcntum, is related to its component transverse 
to the photon beam axis, \k\, as follows: 

q''^P + ip, P = ^2 0^2{e- 1). (10) 

The final-state phase volume may be expressed as 

dr = d^kd%6{{K -ki)^-l)5{{k + ki)^-l)5{{p-k)^-l) 
d^/cd^/ci dx 



4s(l - ak){x + (3k){l - x) 



d^kd^k^dxf (3k\ 

Asx{l-x) i^ + "^-^J- ^''^ 
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In terms of these variables, the total differential cross-section may be written 
in the form: 

„3 



da 



with 



''1212 



''1212 



O1234 



CH278 



and 



a 



^1256 — 



•^1212 + 



«1212l -"Sttfc - "S- 



+ a 



2e 2k'' 

"Ctl234 H ; — 7r^l278 



AB 



1-x 



4:X{1 



Rn{B-Af 



X 
Sl 



1212 



(12) 



-4x(l-x)-^, 



B-A 



{x + l3k)\ AB 
+ 4A^2f(3a;-2) 



A^B^ 



ABJ 



8a;(l — x) 



-^1256 \ <^ki<lx<l'k, 

RijP -A) 
AB2 



Rii{B-A) 



AB'^ 

Z2 



+ (1 - x; 



Rii{B-A) 
A^B 



- 4(3x - 2)i?^ - 4p((6a; - 4)|| + (3 - Ax)^ 



I fB-A 
4 



AB AB 
2x{l-x 



j - 2a;(l 



^RiAB- A) ^ ,^ ^ R 
-^) V ' +2x{l-x)-^ 



Ri{B-A) Ri\ 
AB2 52;' 



^^-^^^^^(.- + (l-2.)-.--(l-.)(---)) 
'^-(^-^)-4(l-x)| + 4(l-.)§-4(2-.)|, 



+ 4- 

1 /B 
'4 



AS 
2x(l - 



AB 

A P 
+ 



X 



AB 

'AB 



+ 2x{l - x) 



Rii{B-A) 
A^B 



(13) 



1 P 

Rn^ekie*ki, Ri ^ -(eki e*k + ek e*ki), R ^ ek e*k ^ —(1 + P). 

In general, the limits of variation for the parameters of the created pair are 
imposed by experimental cuts together with the following relations 



{K — ki)o = w{l — x) > m, 
Sl < s, 



{ki + k)o = {x + l3k)oJ > m, 
s — 2uj 
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or 



2m^ 



<{x + pk,l-x), 0<kl<s{x + /3fc)(l -x) = A. 



3 The inclusive distribution of the recoil elec- 
tron 



The leading and non-leading contributions to the inclusive cross-section in 
recoil-electron momentum may be organized as 



27T 



da 



a 

ro = — , 
m 



where 1(212 comes from the first term in Eq. (0) and corresponds to 
the second. The leading and associated non-leading contributions may be 
decomposed as 



l-e 



'1212 



dx 



d^fci 
'^1212 



(14) 



fc2<A 



l-e 



dx 



d'^h 
a 



TT 



1212 ~ -'1212 ' -'1212- 



0<kl<oo A<fc2<oo 



Using the table of integrals provided in the Appendix, for the first term in 
brackets we obtain 



-'1212 



1 l-e ^2 

y" dz y" dxj h 8x(l - x) 

1 

-_ |d^|_(i_2e + /3fc) + - 



47 

4 + P 
47 



R 



7 



7 



(15) 



where 7 = 1 + 2;(1 — z)/c^. For e = 0, we reproduce the result given in Eq. (^. 
In order to see this, one may use the expansion 



— = 1 - -A;^ H A;^ + . . . 

7 6 30 



and the relation between k"^ and given above in Eq. (|1^). The second 
term, 1^2125 be calculated using the expansion of 0^212 for kf 1 (see 
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Eq. (pOD in the Appendix) 



Tt2 
-'1212 



2k^ 



(L-ln2-l). 



(16) 



The quantity /" may also be expressed as a sum: J" 



'1212 



Consider now the contributions arising from corrections to the leading term 
(see Eq. (|l|)): 



l-e 



■1212 



dx 



TT 



'^1212 ~(^k — 



/3k 



X 



k1<A 



1-t 



dx 



d'h (1 -^)fc 

'^1212 



TT 



1-e 



dx 



x(l — x) 



SX 
^2 



TT 



''1212 



1 + (fci + (l-x)A?)' 



A:2<A 



— _|_ 7'C2 _|_ rc3 

~ -'1212 -'1212 -'1212 • 

The first term on the RHS of Eq. ( [T7| ) gives 



(17) 



-'1212 



^/d.{|(L-ln2-l) + ^ 



4 + fc2 ^z(l-2) 

^ ^ 

47 7 



It is convenient also to present the second term as a sum of two parts. 
The first, containing L^, comes from the 05*212 term, no n- vanishing for both 
X — > and x — * 1: 



TC2 
-'1212 



1-e 



dx 



x(l — x) 



d^A;! xA + {l - x)B - x(l - x)P 



TT 



AB 



dz- 







7 



;(L2-ln2 2) 



TT 



The remaining terms are 

l-e 



TC3 

-'1212 



4 

-- / dx 



A;f<A TT 



xA + (1 - x)5 - x(l - x)fc^ 

























BV 





ri + r2. 



The last term in the first brackets, which gives rise to r2, is ultraviolet con- 
vergent upon integration over ki, 

4fc2 J- r 4 + P z(l -z)l 
r2 = -— dz -1 + —— + ^ Lr , (19) 

whereas the other is 



The non-leading term, 0^2125 along with the interference contribution (see 
Eq. (p!2D ) to the inclusive cross-section yields 

J». = i! /dJ-^(L-l„2)-4 + i±^ + iiil^4 (21) 

s J I 7 7 7 J 

The contribution of the structure ai256 is exactly zero to order 1/s. To make 
this clear, one can bring it into the following form, 

which is antisymmetric with respect to interchange of x and 1 — x. 



4 Conclusions 

Our result for the order 1/s correction to the inclusive cross-section is the sum 
of the expressions for /f2i2^; -^Si? '^^ ^^^^1 /" given in Eqs. (|T^-^T]). Organizing 
the inclusive cross-section in the form 



27r- 



da 



2q 



dq ( 



2(e-l) ^/ 2(£-l)^ 
ao + ^ -ai -P{bo + — -hi 



where gq and 60 are given above in Eq. (0) and 



a-i = - 
2 



-L^ + C- 2LJe + 1) 



C = In^ 2 + 2 In 2 + 



IT 



4 ^ 1.1566, 
9 



\\i{e + q) 




a) 



b) 



Figure 2: The dependence of the correction, A^, to the asymmetry as a 
function of a) s and b) q. 



we extract the asymmetry, 

A = A + AA = — + „ ■ (22) 

ao s ag 

The expansion of can be recast in the form (and this is our final 
result). 



sal 



L'-C-l + Lg{5 ~L' + C)~ 2Li{e + 1) 



(23) 



The the dependence of this quantity on q at fixed s and vice versa is shown 
in Figs. ^ and b (recall that we have set m = 1). 

For small enough s < 10 the terms of order of 1/s", for n > 2, become 
essential and the approach presented in this paper is not applicable. 

We should like to point out that our results are in qualitative agreement 
with those obtained using the HELAS code and reported in the talk 
delivered at the workshop 0]. 
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Appendix 

Since the amphtudes of the gauge-invariant sets of diagrams (3,4), (5,6) and 
(7,8) are suppressed by at least one power of fc^/s as compared to amplitude 
(1,2), we need only consider interference terms. Thus, for the modulus of 
the matrix element, squared and summed over fermion spin states, we have 
within leading (~ s^) and non-leading (~ s) accuracy, in order 

V|,^|2_ (4^^)' f 2Ti234 2Ti256 2ri278 T^^ujl ' ak? ] . . 

^' ' P \ S{l-X) Si SX k2 /' ^ ^ 

where 

^1212 = Tr{(#+ 1)7^(^+1)7,} 

X Tr {(^- ^1 + 1)0^2^(^1+ ^ - l)Oi:} e,e:, 
Ti234 = Tr{(^+1)7;.(^+1)7.(^- ^i + l)0r2'(^i+ ^-l)03:}e,e:, 
Ti256 = Tr{(^+ 1)7^.(^+1)056} (25) 

X Tr {(^- ^1 + l)0r2'(^i+ ^ - 1)7.} e.el, 
Ti278 = Tr{(^+1)7;.(^+ 1)0^3(^-^1 + 1)0^2^(^1+^- 1)7.} e^e:, 

0^2 = -^^7.(^-^i-^+1)7a-^7a(- ^1 + 1)7., 

03^4 = -|7m(^- ^1- ^+1)7a--7a(#-^- ^ + 1)7m, 
n s 

O^J = ^ [-7a(#- ^- ^ + 1)7m + 7m(^+ ^ + 1)7a] , 
Oi;^ = -^7a(- ^i + 1)7m + ^7,.(^+^ + 1)7a, 
and q = p — k. 

The quantities Tjjfe; are related to the aijki given in Eq. (jT^) by 

T1212 = 16sa;(l - x) [a?2i2-5 + a}2i2]> 

2^1234,1278 = 16s^a;(l — a;) ai234,l278) 

^1256 = 16SX(1 - X) ai256- (26) 

To perform the integration over ki, we introduce an ultraviolet cut-off kf < A, 
which may be omitted when calculating convergent integrals for the correc- 
tions. The integrals containing A ot B are (hereinafter we omit the terms of 
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order of 

d% 1 



d^h 1 

fi2 



vr 



52 



vr 



fd'kiRu I,, . f 



d^fcii? 
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TT 



52 



-(lnA-l) + i?. 



(27) 



In order to avoid linearly divergent integrals, we combine denominators as 
follows: 

1 } dz 



AB 







d^kl 


1 


TT 


AB 


d% 


Ri 


TT 


AB 


d% 


Ru 


TT 


AB 



dz 1 , , , 

— = -ln(e + g) = L^, 

7 Q 



zdz 



-R 



7 



' 



(2^ 



dz 



1 1 z^ - 

-(InA- 1)- -ln7 + — 

2^ ' 2 ' 



To evaluate the quantity ri, we use the following set of integrals: 
^\xA+(l-x)B]^ = i?flnA + xP), 

TT B^ ^ ^ 

'^^^'[xA+{l-x)B]R^^^^ = -/?(lnA-l + xP 



TT 



d^ki 



TT 



[xA + (1 - x)B]Rn 



AB^ 
A-B 
AB 



(29) 



«4 



InA 

2 



In deriving ([l6|) , for 0^212 averaged over angles with /c^ 1, we make use of 



''1212lfc2;^l 



(1 -2x(l -x)) 



(A:2)2 

and to perform the angular averaging in ri we use 

Rnikhy likiyf^ + 2R]. 



(30) 
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